Pinning of vortices in a Bose-Einstein condensate by an optical lattice 
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We consider the ground state of vortices in a Bose-Einstein condensate. We show that turning on 
a weak optical periodic potential leads to a transition from the triangular Abrikosov vortex lattice 
to phases where the vortices are pinned by the optical potential. We discuss the phase diagram of 
the system for a two-dimensional optical periodic potential with one vortex per optical lattice cell. 
We also discuss the influence of a one-dimensional optical periodic potential on the vortex ground 
state. The latter situation has no analogue in other condensed-matter systems. 
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Introduction — The effects of a periodic array of pin- 
ning centers on vortices in superconducting materials 
have attracted a lot of experimental 0, Q, IE 0, 0] and 
theoretical 0lllEllE3,IllElEi attention. Of par- 
ticular interest is the effect of the pinning potential on the 
melting of a vortex lattice. The vortex lattice is known 
to melt via a first-order transition in clean materials i 
whereas the presence of pinning centers significantly en- 
riches the phase diagram, due to the intricate interplay 
between the vortex- vortex interactions, pinning poten- 
tial, and thermal and quantum fluctuations 0,0, Si- At 
zero temperature and for strong pinning, the system has, 
depending on the number of vortices per pinning center, 
i.e., the filling factor, various phases where the vortices 
order in a periodic array |9, mi LLll LL2( • If the pinning po- 
tential is weakened, the pinned vortex lattice undergoes a 
first-order transition to a deformed triangular Abrikosov 
lattice [iH- 

Recently, the experimental study of vortices in super- 
fluids has been complemented by the experiments with 
rotating atomic Bose gases 0, 0, 0, . Within this 
field it has become possible to experimentally study the 
dynamics of a single vortex line in great detail |20( , 
leading to an enhancement of the theoretical interest in 
the dynamics of a single vortex in a Bose-Einstein con- 
densed atomic gas [2l| . 

Another interesting development in the field of atomic 
gases is the possibility to trap atoms in a periodic po- 
tential using a so-called optical lattice. Here, one uses 
the dipole force which the atoms experience in an off- 
resonant light field. Using an optical lattice Greiner et 
al. \22\ were able to experimentally observe the transi- 
tion from a superfluid, where the atoms are delocalized 
across the lattice, to a Mott-insulating state where the 
atoms are localized onsite |23j . 

A common feature of the experiments with ultracold 
atomic gases is that these systems are very clean. There- 



fore, vortices in a Bose-Einstein condensate do not ex- 
perience an intrinsic pinning potential and the observed 
vortex lattices are triangular Abrikosov lattices. In this 
Letter, we study the ground state of vortices in a Bose- 
Einstein condensate. We show that turning on an optical 
periodic potential leads to a transition from the triangu- 
lar Abrikosov lattice to phases where the vortices are 
pinned by the optical potential. We restrict ourselves to 
the case with one vortex per unit cell of the two dimen- 
sional optical lattice. 

Interestingly, the precise knowledge of the optical lat- 
tice potential allows for a microscopic and quantitative 
calculation of the phase diagram, as opposed to supercon- 
ducting materials, where the pinning potential is known 
only phenomenologically. In the experiments with ro- 
tating Bose-Einstein condensates, the vortex lattices are 
relatively easy to observe, which allows for a detailed ex- 
perimental study of the transitions between the various 
phases of the vortex lattice as one tunes the strength 
of the optical potential. Moreover, applying a one- 
dimensional optical lattice leads to pinning "valleys" in- 
stead of pinning centers, and, therefore, to pinned phases 
of the vortex lattice which have no analogue in other 
condensed-matter systems. 

Two-dimensional optical lattice — Our starting point 
is the hamiltonian functional for the macroscopic con- 
densate wave function \l/(x), given by 



+ Vol(x)-/x 



— — + -g * x 2 

2m 2 

*(x) . 
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Here, m denotes the mass of one atom which interacts 
with the other atoms via a two-body contact interac- 
tion of strength g = Aira s h 2 /m, with a s the s-wave scat- 
tering length. The optical lattice potential is given by 
Vol(x) = sE^sin 2 (qx) + sin 2 (qy)] with Er the recoil 
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energy, q the wavenumber of the lattice, and s > a 
dimensionless number indicating the strength of the op- 
tical lattice. Note that considering one vortex per unit 
cell of the optical lattice implies that we take the two- 
dimensional vortex density equal to n v — q 2 /Tr 2 . The 
chemical potential that fixes the number of atoms in the 
condensate is represented by /i. 

In the following we consider for simplicity a conden- 
sate with infinite extent in the x-y-plane and tightly con- 
fined in the z-dircction by an harmonic trap with fre- 
quency u> z . The latter assumption allows us to neglect 
the curvature of the vortex lines, and leads effectively 
to a condensate thickness d = \J irh/(muj z ) in the z- 
direction. In the Thomas-Fermi limit, where we neglect 
the kinetic energy of the condensate atoms with respect 
to their mean-field interaction energy, the global density 
profile of the condensate in the optical lattice is given by 
71tf(x) = n - [Vol(x) - sEn]/g, with n = \p - sE R ]/g 
the average density of the condensate. To find the poten- 
tial energy of a vortex in a Bose-Einstcin condensate in 
an optical lattice as function of its coordinates (u x ,u y ), 
we use the variational ansatz for the condensate wave 
function 



*(x) = Vn TP (x) 0[|x - u|/f - 1] exp[#(x, u)] , (2) 



with £ = l/^/87ra s n the healing length that sets the size 
of the vortex core, 0(x, u) = arctan[(y — u y )/(x — u x )\ 
the phase of the vortex, and Q(z) the unit step function. 
For the above ansatz to be valid, we have assumed that 
the vortex core is much smaller then an optical lattice 
period, q£ <C 1, and that the strength of the potential is 
sufficiently weak, sEr <C /i. Note that the above varia- 
tional ansatz indeed describes a vortex along the z-axis 
at position (u Xl u y ) in the x-y-planc. 

Substitution of the ansatz in Eq. (J2J) in the hamiltonian 
functional in Eq. and isolating the contribution due 
to the presence of the vortex leads to [24| 

C^2d(u) = -^—sE R Q(q£)[cos(2qu x ) + cos(2qu y )} . (3) 
oa s 

Here, we defined Q(z) = J 1 (2z)/(2z) + dpJ (2zp)/ p, 
with Jq and J\ the zeroth and first order Bessel function 
of the first kind. It is clearly seen that the potential en- 
ergy is minimal if the vortices are located at the maxima 
of the optical potential. This is expected, since at these 
maxima the condensate density, and hence the kinetic en- 
ergy associated with the superfluid motion, is minimal. 
The expression in Eq. (HJ is the pinning potential expe- 
rienced by vortices in a condensate loaded in a optical 
lattice. If the pinning potential is sufficiently strong and 
we have one vortex per unit cell of the optical lattice, the 
ground state is a configuration in which each vortex is 
trapped or pinned by an optical lattice maximum. For 
a two-dimensional optical periodic potential, this is the 
square pinned lattice (SP), shown in Fig. ^ 



To determine the phase diagram in detail, we calculate 
the total energy per vortex for an arbitrary vortex lat- 
tice. This approach neglects the fact that for very weak 
pinning potentials the triangular Abrikosov lattice will 
be slightly deformed . Generally, the vortex lattice is 
parametrized as follows 



7T / VI + a/3 
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(4) 



with ni S Z and < a, (3 < ^. The transformation ma- 
trix conserves the area of the unit cell, and thus ensures 
that we are considering configurations with equal vortex 
density. The more familiar parameters of a unit cell of a 
two-dimensional lattice, the angle ip and the ratio of the 
length of the sides of the cell n = L1/L2, are related to 
a and /3 by 



cosp (a + /3)V1 + a P sinp 
k 1 + a(3 + a 2 k 



1 



1 + a/3 + a 2 



(5) 



The pinning energy per vortex is found by substituting 
Eq. I0J in Eq. iJSJ, summing over all n^, and dividing the 
result by the number of unit cells. In the limit ni — > 00 
we find 



E&(a,0) = -—sE n Q(q£)[6p t0 + S afi ] 



(6) 



Our next task is to determine the interaction energy 
per vortex. In two dimensions, singly-quantized vortices 
with equal orientation experience a logarithmic long- 
range interaction V(r) = —2Tidh 2 n/m\og(r/£) for r ^> £ 
|25|. The interaction energy E- lnt per vortex for an infi- 
nite two-dimensional lattice of vortices was calculated by 
Campbell et al. 26}. Cast in a dimensionless form, their 
result reads 



Et 



E„ 



7T Sllip 
6 K 



(Trh 2 dn/m) 
log{n~ 1 [l-2e- 27rj|sinvl/K 

e -4iu"|sint/>|/K j j 



log 



2tt 



sin p 



(2nj 



cos p 



(7) 



It is important to realize that the interaction energy per 
vortex is divergent for an infinite lattice, and that the 
above expression gives the relative interaction energy for 
configurations with equal vortex density. 

The expression for E- mt together with the expression 
for the pinning energy in Eq. |0 enables us to mini- 
mize the total energy as a function of a and (3, and to 
determine the vortex lattice ground state for different 
strengths of the periodical optical potential. Performing 
this procedure as function of the dimensionless variables 
ql; and sEr//i leads to three phases, which are physically 
due to the the competition between pinning and interac- 
tions. For a very weak optical potential we find that vor- 
tices arrange themselves on a triangular Abrikosov lat- 



tice (AB), with a = f3 = Jl/\/Z— 1/2, as expected. In 
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FIG. 1: Vortex phase diagram of a Bose-Einstein condensate 
in a two-dimensional optical square lattice, with one vortex 
per unit cell of the optical lattice. Three phases are relevant: 
a square and fully pinned configuration (SP), a triangular 
configuration where half of the vortices are located at the 
pinning centers (HP), and the unpinned triangular Abrikosov 
vortex lattice (AB). 
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FIG. 2: Vortex phase diagram of a Bose-Einstein condensate 
in a one-dimensional optical lattice. Two phases are relevant: 
a pinned triangular configuration (PT), and the unpinned 
triangular Abrikosov vortex lattice (AB). We calculated the 
phase boundaries for £ = .Oln/qo (+), £ = . 005^/^0 (x) and 
£ = .00l7r/qo (*)• Note that at the line q/qo = 1 the phases 
coincide. 



this phase, with E- lnt = —1.32112, the interactions domi- 
nate over the pinning. On the other hand, when the pin- 
ning energy dominates over the interaction, we find the 
square pinned lattice (a = j3 = 0) 0. This phase has 
E[ nt = —1.31053. In the intermediate regime, where the 
pinning and interactions are equally important, we find a 
phase in which half of the vortices are pinned (HP) [l3| . 
and the lattice has a triangular structure (a = | , P = 
and Ei nt = —1.31800). In the zero-temperature phase di- 
agram, shown in Fig.^ the different phases are separated 
by first-order phase transitions and the phase boundaries 
are given by 

(sE K \ _ .00623 fsE R \ _ .01494 

V f 1 /AB-HP A 4 / HP-SP 

(8) 

One- dimensional optical lattice For a one- 

dimensional optical lattice the single vortex pinning 
potential equals J7id(ii) = -^-sE R Q(q£ i ) cos(2qu x ) 
and hence the minima of U±r> act as pinning "val- 
leys" . The pinning energy per vortex reads E^ = 

— ■^-sEu,Q(q£,)Sp_o within the parametrization in Eq. 
Two phases are distinguished in this case, i.e., a pinned 
triangular lattice (PT), shown in Fig. |3 and the un- 
pinned Abrikosov vortex lattice. Interestingly, the inter- 
actions always favor a triangular lattice since the vortices 
are allowed to arrange freely in the y-direction. 

Consider now the case that the wavenumber of the op- 
tical lattice is such that the AB lattice and the PT lattice 
coincide, i.e., qo = 2ir^/3/ (3L), with L the intervortex 
distance. If this configuration is disturbed by changing 
the optical lattice wavenumber to arbitrary q there will 
be a competition between the AB lattice and the PT lat- 
tice. The unit cell of the latter is, for wavenumber q and 



at equal vortex density, described by 




The interaction energy per unit cell is found by substi- 
tution of Eq. © in Eq. (JJJ) . From this we find the zero- 
temperature phase diagram, depicted in Fig-EJfor various 
values of the healing length. The generic behavior is such 
that for given strength of the optical lattice and for small 
deviations from q , the vortex lattice stays pinned, i.e., 
the vortices are dragged along with the pinning valleys. 
At certain q the phase boundary is crossed. Then the 
vortices "jump" back to their original positions, forming 
an Abrikosov lattice again. 

Discussion and conclusions — Since we have consid- 
ered the ground state of vortices, we have implicitly as- 
sumed that the optical lattice is co-rotating with the 
Bose-Einstein condensate. Although this is very difficult 
to achieve experimentally, it has, however, recently been 
proposed to create effective magnetic fields, and there- 
fore effective rotation, by optical methods [23, |28| • We 
believe that the calculations presented here are relevant 
for such a situation. 

Of particular experimental interest are the collective 
modes supported by the pinned vortex lattices. In the 
absence of an optical lattice potential, the dispersion re- 
lation for the gapless Tkachenko modes has been mea- 
sured by Coddington et al. 29] . We expect that these 
modes acquire a gap in the presence of a periodic optical 
potential due to the fact that the translational symmetry 
is not broken spontaneously, but by the optical lattice. 
In the case of the SP lattice this gap is easy to calculate, 
since the zero-momentum Tkachenko mode corresponds 
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in this case to a simultaneous in-phase precession of all 
the vortices around the maxima of the optical lattice. 
Hence the gap is given by hui p , with uj p the precession 
frequency |2l| 



irhq 2 Q{qOsE R 



mgn 



(10) 



In future work we intend to study the collective modes 
of the pinned vortex lattices in great detail. 

In this Letter we have restriced ourselves to the case 
of one vortex per unit cell of the optical lattice. Another 
direction for future work will be a detailed study of the 
pinned phases at other filling factors, where other types of 
pinned vortex lattices are known to occur 0, 0, Oj], 0] . 
Since we have only been considering infinite lattices, we 
intend to study also the effects of the finite size of the 
system, which may be significant |30j|. With respect to 
this latter remark it is also important to note that the 
harmonic magnetic trap used in the experiments induces 
an additional feature in the density profile of the conden- 
sate that may have important effects j^H, and therefore 
also requires further study. 

Apart from these interesting possibilities, yet another 
direction would be to consider more strongly-correlated 
regimes that occur at fast rotation, and to study the ef- 
fects of the periodic optical potential on the melting of 
the vortex lattice [j^]. One would expect that in this 
regime the effect of quantum fluctuations, i.e., quantum 
tunneling of the vortices through the potential barriers 
of the pinning potential, becomes important. Conversely, 
in the limit of a strong optical periodic potential, an in- 
teresting topic is to study the effects of the rotation on 
the Mott-Insulator transition In conclusion, rotat- 
ing Bose gases in an optical lattice provide an interesting 
system to study new quantum phases of matter jsij, as 
well as phenomena known from condensed matter in a 
new context. 
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